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I.  INTRODUCTION  ANl>  SUMMARY 

Sport  t  all  v  r  esolved  infrared  r  adiance  profiles  of  Ihe  Earth's  limb  can  be  in¬ 
verted  to  obtain  vertical  distr  ibutions  of  temperature  and  species  concentration. 

The  atmosphere  below  approximately  70  km,  where  l.TE  prevails  and  the  undis¬ 
sociate  I  CO.,  distribution  is  kne;  ,  can  be  probed  using  broadband  limb  measure¬ 
ments.  as  demonstrated  bv  House  and  Ohring,  1  Cl i lie*  et  al,  ^  Russell  et  al,  ^  and 
others.  The  added  dimension  of  spectral  resolution  offers  the  potential  of  obtaining 
unambiguous  information  in  the  non-l.TE  altitude  regime. 

The  limb  spect  ral  radiance  profile  cannot  in  general  be  expressed  as  a  linear 
'unction  of  the  vertical  temperature  and  concentration  profiles.  Consequently,  the 
inversion  must  generally  be  performed  bv  some  sort  of  recursive  algorithm. 

/itch  >r  et  al  list'll  a  recursive  direct  nonlinear  technique  to  recover  vertical  pro¬ 
xies  if  excited  nitric  nxi  !«•  concentration  and  kinetic  temperature  above  a  100  -  km 
altitude  Iron  spectral  data  measured  bv  the  rocket-borne  SHIRK  C\  K  spectrum  - 
a  e  i  ,  These  inversions  used  SHIRK  lata  tor  i  lie  \()  fundamental  hand  near 

fi.  (  im  ii.i  i -able  results  i  mid  have  been  ihtained  bv  m  approximate  linear 

i  ..  Hi  i  I  In  au.-e  j round  state  \<  >  is  optically  thin  for  tangent  paths  above  -  100  km 

(Re,  (  I  '■  I  '  II'  |  III  111  I  C  It  I  O' |  Se  111  be  r  I  Ui'ii!  I 

i  I  me  u>  :  ■  i"  !  a  r  ,'i  mm  her  i :  re  le  rein  es  ■  i  ted  a  ti  ne,  t  lie  v  will  not  he  listed  here. 

S,,.  |{.  te-en  e  Il.il>.  |;>.  » 


'  I 


for  the  fundamental  band.  That  is.  the  measured  radiance  profile  and  the  excited 

state  density  profile  have  a  relationship  that  is  approximately  linear  for  this  special 

case.  A  nonlinear  method  can  be  adapted  more  easily  to  the  general,  non-optically- 

4 

thin  case,  which  is  why  we  chose  to  use  the  nonlinear  method  in  our  earlier  study, 
referred  hereafter  to  as  ZSNS. 

The  direct  nonlinear  method  starts  with  an  initial  guess  for  the  solution  profile, 
for  example,  the  excited  NO  vertical  density  distribution.  The  corresponding  limb 
radiance  profile  is  computed  and  then  the  guess  is  relaxed  according  to  some 
measure  of  the  disagreement  between  the  computed  and  measured  limb  radiance 
profiles.  The  relaxation  equation  used  in  the  earlier  study,  Rq.  (8)  of  ZSNS,  re¬ 
vises  the  guess  according  to  ratios  of  measured  and  computed  radiance  values. 

The  process  is  iterated  until  the  solution  profile  has  converged. 

The  purpose  of  this  studv  is  to  examine  further  the  inversion  capabilities  of  the 
relaxation  method  used  in  ZSNS,  and  to  determine  its  sensitivity  to  noise.  We  were 
interested  particularly  in  answering  the  following  questions: 

(a)  How  accurately  can  the  method  recover  sharplv -peaked 
distributions?  In  the  limb  viewing  geometry,  a  laver 
far  below  the  peak  is  observed  through  much  denser 
lavers  above  it. 

(b)  How  accurately  can  the  temperature  profile  be  recovered 
when  there  is  noise  in  the  data''  The  temperature 

profile  recovered  in  ZSNS  had  large,  noise -like  excursions. 

These  questions  were  answered  bv  generating  synthetic  spectral  limb  radiance 
profiles  and  then  inverting  them,  both  before  and  after  the  addition  of  normally 
distributed  pseudo-random  numbers  representing  noise.  We  used  the  NO  funda¬ 
mental  (Av  -  1)  band  and  assumed  that  ground  state  NO  is  optically  thin  for  the 
complete  range  of  tangent  paths  involved  in  the  simulations,  which  corresponded 
to  tangent  heights  between  85  and  130  km.  Thus,  the  equations  used  for  the  direct 
computation  of  radiance,  spectral  radiance  and  the  relaxed  guess  are  those  given 
in  ZSNS.  The  inversion  method  modified  to  handle  the  non -opti call v -thin  case  (and 
to  recover  ground  state  as  well  as  excited  state  density)  would  be  expected  to  have 
a  slightly  higher  sensitivity  to  noise  than  will  be  indicated  bv  the  present  results. 

In  the  noise-free  simulations,  the  excited  NO  vertical  density  profile  was 
represented  bv  a  Gaussian  function  and  also  bv  a  double -peaked  distribution  which 
is  the  sum  of  two  Gaussians.  The  simulations  with  noise  used  onlv  the  double -peaked 
distribution.  The  added  white  spectral  noise  bad  two  different  tons  levels,  corres¬ 
ponding  to  maximum  S  Vs  m  500  and  100;  these  arc  the  signal  -lo-noise  \  a  lues  at 
'he  .  ombination  if  wavelength  anl  tangent  height  that  pro  duces  the  maximum  limb 
spectral  radiance.  Idle  -oectral  resolution  was  the  same  as  the  SI’li.’t'  instrument, 
approximately  0.05  pm,  who  It  resolves  the  \( )  band  into  13  soectral  <■!<  moots. 


♦ 

The  excited  NO  density  profile  NO  (H)  was  retrieved  by  inverting  the  synthetic 
limb  profile  of  band  radiance  (integrated  spectral  radiance),  whereas  the  tempera¬ 
ture  inversion  used  13  profiles  of  spectral  radiance.  As  in  the  ZSNS  study,  the 
NO  (H)  solution  and  the  set  of  solutions  obtained  by  inverting  each  spectral  radiance 
profile  are  combined  to  obtain  the  profile  R(H.  A),  the  normalized  spectral  emission 
per  unit  volume  as  a  function  of  altitude.  The  normalized  spectral  emission  R(T,  A) 
was  also  computed  for  various  temperatures  T.  A  systematic  comparison  of  the 
retrieved  Kdl,  A)  and  computed  K(T,  A)s  yields  the  temperature  profile  T(H).  The 
initial  guess  in  each  of  the  inversions  consisted  of  a  constant  (al titude -independent ) 
a  l  lie. 

1  he  results  of  the  simulations  are  displayed  in  graphs  in  the  following  section. 
Kxaminali  >n  of 'he  results  leads  to  the  following  conclusions: 

(ai  In  the  noise-free  simulations  the  retrieve*!  NO  (III  and  1(11) 
were  ucar  iv  identical  to  the  models  used  to  generate  the 
1 1 1 1 1  *  ■  1 1  •  ■  limb  radian  ornfilcs,  e\i  cpi  at  the  verv  lowest 
i !  * :  i  it where  the  *  -  x  *  *  i  t « ■  *  I  i«  *n  i  t  v  w  as  two  >rders  of 
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(e)  Examination  of  the  results  indicates  that  the  random  peak 

excursions  in  the  retrieved  temperature  profile  at  the  altitudes 
of  the  NO  peaks  are  approximately  ±5°C  for  maximum  spectral 
S/N  =  500  and  approximately  ±25°C  for  (S./N.)  =  100.  The 

temperature  profile  that  we  retrieved  from  the  SPIKE  NO  data 
(see  ZSNS)  has  excursions  of  ±7  5°C  near  the  altitude  of 
maximum  NO  .  approximately  130  km.  This  would  be  consistent 
with  the  present  simulation  results  if  the  spectral  signal -to -noise 
of  the  SPIRE  data  for  tangent  height  =  130  km  was  approxi¬ 
mately  30  (after  the  data  is  averaged  over  0.05  pm  intervals, 

as  done  by  ZSNS).  Examination  of  the  SPIRE  data  for  this 
tangent  height  shows  that  the  enhanced  S^/N^  is  somewhat 
greater  than  30,  but  probably  not  more  than  60.  Hence,  the 
present  results  support  the  conclusions  by  ZSNS  that  detector 
noise  was  the  major  cause  of  the  noise-like  excursions  in  the 
temperature  profile  retrieved  from  the  SPIRE  NO  data. 

(f)  The  results  emphasize  the  importance  of  achieving  high  signal - 
to-noise  in  limb  spectral  radiance  measurements,  and  also 

of  considering  the  tradeoffs  between  spectral  resolution  and 
signal -to-noise.  This  study  did  not  consider  the  effects  of 
smoothness  constraints  on  the  inversion  solutions,  which  might 
reduce  considerably  the  S/N  requirement. 


2.  RESULTS  OF  SIMULATION 


The  single -peaked  Gaussian  distribution  used  to  represent  the  excited  density 
profile  NO  (H)  in  the  simulations  is  given  by 

NO*(H)  -  5  X  104  exp  [ ( II -  1 1 5 ) / 5 ] 2  ,  (1) 


+  -:t 

where  NO  has  the  units  cm  and  (lie  altitude  H  is  in  kilometers.  The  density 
profile  peaks  at  115  km  and  is  e  *  times  the  peak  value  at  110  and  120  km  altitude. 
The  double -peaked  distribution,  given  bv 


(2) 


has  a  relative  minimum  at  lOi.a  km,  where  NO  is  approximate! v  one-tenth  the 
maximum  value  (-5  ■  10*  cm  >  which  on  urs  at  100  km  and  115  km. 
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The  translational  temperature  profile  used  in  the  computation  of  synthetic 

7 

radiances  is  the  profile  given  by  the  Jacehia  model  for  an  exospheric  temperature 
of  1050K  (based  on  the  U.  S.  Standard  Atmosphere  197  6).  The  profile  is  displayed 
in  graphs  that  compare  the  model  T(H)  to  the  inversion  results,  which  are  des- 

JjC 

cribed  below.  Figures  1  and  2  show  the  two  different  NO  (H)  models  and  the 
corresponding  computed  limb  profiles  of  band  radiance.  The  13  computed  limb 
spectral  radiance  profiles  are  not  shown. 

Figure  3  compares  the  single -peaked  "actual"  or  model  NO  distribution  and 
the  cor  responding  retrieved  profile  when  no  noise  has  been  added  to  the  synthetic 
band  radiance  profile.  However,  the  synthetic  band  radiances  were  rounded  to 
four  significant  figures  before  they  were  inverted.  Note  that  the  error  is  sig- 

if 

nificant  only  at  the  lowest  altitudes,  where  NO  (H)  is  less  than  0.01  times  the  peak 

if 

NO  .  Figure  4  shows  that  the  error  has  essentially  the  same  behavior  for  the 
double-peaked  case. 


RADIANCE  IWcm^sr'1) 


N0*!H)  tm"3t 

if 

Figure  I.  Single -peaked  Gaussian  Model  for  NO  (H)  and 
tiie  Limb  Profile  of  Hand  Radiance,  Determined  by  NO*(II) 
and  the  Model  Temperature  Profile  T(H) 


■Picchia,  !..  (1977)  Thermospheric  Temperature,  Densitv  anti  Composition: 
N'’«  Models.  Smithsonian  Astrophvsi cal  Observatory  Special  Report  375. 


Figure  4.  Same  as  Figure  3  Except  that  the  Actual 
.\0*(H)  Is  Represented  by  the  Double -peaked  Model 


Figure  5  compares  the  actual  or  model  T(H)  to  the  temperature  profiles  re - 
rievrd  lor  the  single-  and  double  -  peaked  cases,  when  only  rounding  noise  is 
oresent.  The  errors  in  retrieved  T(H)  are  significant  at  the  same  altitudes  for 
a  hit  h  the  retrieved  NO  (H)  has  significant  errors,  that  is,  when  NO  (H)  is  less 
him  0.01  times  the  (teak  NO  . 

Figure  6  shows  the  retrieved  double -peaked  NO  vertical  profile  when  the 
:eaxitttuc  -  neutral  S/N  is  500.  The  retrieved  profile  is  in  good  agreement  with 
he  i  ord  profile  when  NO  (II)  is  greater  than  '-0,01  to  0.02  times  the  peak  NO 
For  a  maximum  spectral  S/N  of  100  the  actual  and  retrieved  profiles  are  in  good 
igi'cemi'Tit  when  NO  f  1 1 )  is  greater  than  approximately  0.02  times  the  peak  NO  ; 
sec  Figure  7.  However,  at  107.5  km  where  NO  till  has  a  relative  minimum,  the 
relative  error  is  approximately  40  percent. 


Figure  8.  Retrieved 

Temperature  Profile 

for  (S  /N . )  =500 

A  A  max 

Double-peaked  Case 


r>  0  iOi  3«  401' 


Figure  0.  Same  as 

Figure  ii,  Fxeept 

(S  'N  )  -  100 

A  A  m. in 
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:t.  Disci  ssion 


It  ,rs  tn.it  tin'  nonlinear  inversion  technique  investigated  in  this  study  is 

no  less  .ivvurate  below  the  peakls)  in  concentration  distribution  than  it  is  above  the 
pe.i.ki.'),  . i it ; i« > lilt h  the  t  rrors  tend  to  be  more  systematic  at  the  lower  altitudes 
be!ow  a  single  or  double  peak. 

Hosed  on  tin-  results  of  this  study,  it  is  not  surprising  that  the  temperature 
profile  retrieved  from  the  SiTKF  NO  spectral  data  by  ZSNS  show  noise-like  excur¬ 
sions  of  i  75  °C.  The  maximum  S./N  was  between  50  and  ti  0  for  the  .SPIKE  data. 

A  A 

The  present  study  implies  that  maximum  spectral  S/Ns  of  the  order  of  500  are 
required  to  ensure  acceptable  temperature  solutions. 

Translational  temperature  profiles  in  the  thermosphere  have  a  characteristic 
smoothness  exemplified  by  the  model  Till)  used  in  the  simulations.  Based  on  the 
a  priori  knowledge  (or  assumption)  of  "smoothness",  the  retrieved  temperature 
profile  shown  in  Figure  0  would  be  judged  "unacceptable".  So-called  smoothness 
constraints  can  be  incorporated  into  the  inversion  algorithm  to  significantly  reduce 
the  S.  \  requirement.  For  the  double-peaked  NO  cases  modelled  in  this  study,  it 
is  estimated  that  the  S/N  requirement  can  be  reduced  by  a  factor  as  large  as  two 

or  t  h  ree. 

Fee  st ud v  results  indicate,  nonetheless,  the  importance  of  considering  the 
-van  ible  t,.  , deoils  between  spectral  resolution  and  signal-to-noise.  High  spectral 
resolution  might  he  required,  for  example,  when  two  or  more  species  have  ovor- 
!  ppmg  spectral  b  inds,  or  warn  the  available  data  interpretation  algorithm  require 
i  -  i-l  a  ion  of  optic  ll\ -tinn  por*  ions  of  a  strongh  emitting  band.  \\  lien  those  con- 
i . : ■  :■  t  ion-,  do  not  ppK  one  should  opt  for  the  minimum  spectral  resolution  that 
.  i  t  I  o  a  li  .  m  i  .  ontoii  r  I*  rails  I .  d  iona  1  ti-mpe  l*at  u  re ).  The  t  radeoff  s  obviously 
la  o  i ;  .  ,.•••  "r.i'lex  .'.ii  u  t ;  .  i  list  rn  1 1 .  etit  must  recover  s  in  i  u  It  a  neo  us  1  y  the  dis- 
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